A Numerical Simulation of Chern-Simons Inflation 



c/3 



Annie Preston 1 , David Garrison 2 and Stephon Alexander 3 

1 Department of Physics, Harverford College, Harverford, PA 19041 

2 Physics Department, University of Houston Clear Lake, Houston, TX 77058 

3 Department of Physics and Astronomy, Dartmouth College, Hanover, NH 03755 

1} ■ Abstract. In this work, we present results of numerical simulations of the Chern- 

Simons Inflation Model proposed by Alexander, Marciano and Spergel. According 
to this model, inflation begins with a fermion condensate interacting with a gauge 
field. Crucial to the success of this mechanism is the assumption that the Chern- 
Simons interaction would drive energy from the initial random spectrum into a narrow 
band of frequencies at superhorizon scales. In this work we numerically confirm this 
expectation. These gauge fields, when combined with the Friedmann equations, were 
broken into a system of hyperbolic equations and numerically simulated with a novel 
,13 ■ relativistic MHD code. We show that the amplification of horizon sized gauge fields 

produces the conditions to cause cosmological inflation and that the onset of inflation 
are robust against certain fine tunings in the initial conditions. 
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1. Introduction and Motivation 

The paradigm of the early universe is based on the phenomenon of cosmic inflation 
(Guth 1981). It is well known by now that inflation resolves most of the problems of the 
standard big-bang scenario, such as the horizon problem. Moreover, inflation predicts 
a nearly scale-invariant power spectrum of density perturbations which is observed in 
the CMB. Although, there many successes of inflation, there are hundreds of scalar field 
models of inflation which have proven difficult to distinguish with data. Moreover, scalar 
field driven inflation are frought with conceptual and technical issues. Recently Chen 
and Wang demonstrated that the primordial power-spectrum in scalar field inflation 
have anomalous sensitivity to high energy physics (Chen & Wang 2012). These issues 
have motivated the authors (Alexander et al. 2011) to provide an inflationary mechanism 
that is driven by an interaction between vector and spinor fields, otherwise known as 
Chern-Simons inflation. 

Alexander, Marciano, and Spergel propose a model in which the early universe is 
dominated by a gauge field that interacts with a fermion current. This interaction results 
in an equation of state with consistently negative pressure, the condition needed for 
inflation (Alexander et al. 2011). In this model, the gauge field begins as a random, white 
noise spectrum; the authors assume that this evolves into a spectrum of superhorizon 
modes. Because of the complexity of the differential equations, they were not able 
to show this evolution analytically. The purpose of this investigation is to determine 
numerically whether the gauge field transforms from a random spectrum into a spectrum 
of low-frequency modes during inflation. 

Another goal of this project is to establish the robustness of the model. One 
criticism of inflation is that specific theories can focus too much on formalism and lack a 
clear connection to physical processes that actually could have occurred (Kinney 2009). 
Therefore, we aim to strengthen the physical interpretation of an inflationary theory 
while assessing its accuracy; we use a modified magnetohydro dynamic (MHD) code 
based on the Cactus framework to simulate the early universe. Though this is perhaps 
the first cosmological MHD inflation simulation, there have been many astrophysical 
MHD simulations, and there has been a wide range of other types of inflation simulations. 
For example, in (Dubois & Teyssier 2008), the authors use an MHD solver to simulate 
the environment within galaxy clusters. Examples of other inflationary models that 
have been simulated include a back-reacting cosmology with inhomogeneous initial 
conditions (Kurki-Suonio et al. 1987), pre-Big Bang inflation (Maharana et al. 1998), 
and an examination of topological defects due to inflation (Kasuya & Kawasaki 1998). 
Our work presents another method of testing new inflationary theories in an accessible 
way, and introduces the concept of treating the early universe as an MHD like system. 

In Section 2, we provide necessary mathematical background, then we explain 
motivations for all possible mechanisms of inflation. In Section 3, we describe the 
analytical and computational methods used to simulate the specific inflationary theory 
from (Alexander et al. 2011). In Section 4, we present the results of the simulation, and 
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in Section 5, we discuss the results and future research directions. 
2. Inflation Models 

The accelerated expansion during inflation requires that the universe be dominated by 
a negative-pressure component. In many inflationary models, a scalar field dominates 
the early universe and provides the necessary negative-pressure equation of state. The 
following discussion is closely based on (Ryden 2003). This scalar field, 0(f, t), is known 
as the inflaton field. It has a potential V(4>) that depends on the temperature of the 
radiation dominating the universe prior to inflation (Montani et al. 2011). The energy 
density and pressure of the inflaton field are 

*=^+v W (i) 

P*=\&-VW (2) 
If the change in the inflaton field over time is very slow: 

2 < Hc 3 V((f)) (3) 

then we have 

€0 « -P4, (4) 

which meets the requirement for negative-pressure domination. We can see here that 
if there is a period in which <fi is small, and V((f>) is large enough to dominate the 
energy density of the universe, then the inflaton field can cause exponential expansion. 
This model is known as slow roll inflation. However, it is not immediately clear what 
these conditions physically correspond to in the early universe. We can determine the 
necessarily small value of <fi by considering the inflaton field's fluid equation: 

e* + 3if(f)(e* + P*) = (5) 

This gives 

e + 3if (*)e* = -ftc 3 ^", (6) 

which is an equation of motion in which the accelerating force is and the "frictional" 
force, proportional to speed, is 3H(t). This force, known as the "Hubble friction," slows 
the change in the inflaton field until its potential, V ((/>), is at a minimum. When the 
acceleration <fi reaches zero, we have 

■ Ke> dV 

^ = -3W (7) 
We know that the requirement in Equation 3 is necessary for the inflaton field to cause 
inflation, so imposing that requirement gives 

,dV. 2 9H 2 V 
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When the potential energy of the inflaton dominates, which is necessary for inflation, 
the Hubble parameter is 

R = ( 8^ )V 2 = (^)V=» (9 ) 
and our requirement becomes 

where Ep is the Planck energy. Note that during inflation, the Hubble parameter is 
constant. We have found that the inflaton's potential must satisfy the equation above, 
and that its energy density must be large enough to dominate that of the universe, in 
order for exponential expansion to occur. These conditions imply a nearly constant 
energy density of the inflaton field: p ~ 2 + V(0); this requires its kinetic energy to be 
insignificant compared to its potential energy. 

To see how inflation solves the Flatness problem, consider yet another version of 
the Friedmann equation: 

If the scale factor, a, increases exponentially while the energy density, p, stays nearly 
constant, Q must approach 1, since the right side of the equation is constant. During 
inflation, the energy density of the universe is dominated by that of the inflaton field, 
which remains constant just as the energy density of dark energy does. Therefore, 
inflation drives the universe towards flatness. Because the Hubble parameter is fixed 
during inflation, the curvatures immediately before and after inflation are related by 

n f - 1 = (fti - l)e- 2N (12) 

where N is the number of e-folds. According to the Standard Cosmological Model, the 
quantity Q - 1 has increased by a factor of about 10 60 since the end of inflation; therefore, 
the current observational constraint on |Q — 1| gives 

\n - 1| = 10 60 |^- l\e~ 2N < 1(T 2 . (13) 

The minimum number of e-folds needed to resolve the flatness paradox is then 

iV>71 + (14) 

We see that even with significant initial deviations from flatness, enough e-folds of 
inflation efficiently flatten the universe. 



2.1. Chern-Simons Inflation 

Because these scalar field models are widely found and difficult to distinguish from one 
another, Alexander, Marciano and Spergel suggest an alternate model in which a gauge 
field interacts with fermions in the early universe to produce an effective scalar field 
that generates inflation (Alexander et al. 2011). 
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In this gauge field model, we have an energy density 

P=- S J- + I*-JI (15) 

— # 

— » • — * — * — * 

where the gauge field has an electric field, E = A, and a magnetic field, B = V x A. 
We can see that this energy density has an electromagnetic component, which scales as 
a -4 , as well as a gauge field-fermion interaction. As shown in (Alexander et al. 2011), 
if the second term dominates over the first term, while remaining nearly constant, it 
is possible for inflation to occur. This will give us an energy density p = \A ■ J\ and 
pressure P = -\A ■ J\ , corresponding to an equation of state w = -1, which is sufficient 
to cause inflation. 

Just as we gained insight through the scalar field's equation of motion for the general 
case, we will consider the equation of motion of the gauge field in this case. Its action 
is: 

f Jl/t2 pi 1 

S = / d 4 x^(^p - l -d,dd»d + £R9 2 - - A Tr[F aP F^} 
Jm 4 8tt 2 4 

+ lk Tr[Fa ' Fa " ] + V Tr ^ J ^ (16) 
By varying this action with respect to the gauge field, we find the gauge field's equation 
of motion. In terms of Fourier modes, this equation of motion is: 

A h + k 2 A h = -hkA h ^- + a 4 J h (17) 

where the fermion currents are 

V2J h = Ji + hij 2 = ($( 7l + hi^hs^/a. (18) 

Alexander et al. show that the fermion currents have a background solution in which 
they scale as Jo /a . Since they are being suppressed by a factor of a, we can define a 
set of constants Jh where Jh = Jh/a (Alexander et al. 2011). The equation of motion 
for the gauge field is now 

A h + k 2 A h = -hkA h — + a 3 J h (19) 

According to the solutions for the gauge field (Alexander et al. 2011), a value of k < 
corresponding to long-wavelength modes, will result in its exponential growth. However, 
the presence of short wavelength modes, k > -jj- , will result in continued linear growth 
due to back-reaction of the gravitational field. 

Because the universe was a plasma before last scattering, we look to plasma physics 
for an analogy. Specifically, we consider magnetohydrodynamics (MHD), the study 
of electrically conducting fluids. When fluids in motion are electrically conducting, 
currents in the fluid produce magnetic fields that induce forces and thus change the 
dynamics of the system when the system becomes turbulent (Shebalin 2010). This kind 
of turbulence is analogous to the dynamical system changes that could be capable of 
generating inflationary behavior. Computer codes that simulate MHD use the equations 
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of fluid dynamics, along with Maxwell's equations, to analyze the dynamics of a plasma 
system. Simulations such as these are commonly known to lead to what is referred to as 
an "inverse energy cascade". In a hydrodynamic turbulent system, as a system evolves 
in time, energy tends to "cascade" from the longest wavelength modes to the shortest 
wavelength modes as turbulent eddies get smaller and smaller. In MHD systems, the 
reverse happens. Energy cascades inversely from high frequency modes to low frequency 
modes. 

Modifying an MHD code by adding the coupled partial differential equations that 
describe the very early universe could allow us to see the resulting exponential increase 
and then "calming down" of the scale factor. Additionally, a simulation of this type 
will allow us to see directly the physical effects of fine-tuning the initial conditions for 
these coupled differential equations. We can input the equations of motion for the fields 
that are interacting to generate inflation, along with the Friedmann equations, which 
describe the expansion of space. The Friedmann equations will allow us to relate the 
evolving pressure and density of the system to the evolution of the scale factor. 



3. Methods and Equations 

For the numerical calculation, we use natural units but later evaluate the data in terms 
of SI units so that the results can be easily compared to the established values. The 
relevant equation from (Alexander et al. 2011) is the equation of motion for the gauge 
field, A: 

A - V 2 A+ -^-(V x A) - a 4 J, = (20) 

Here, J is the fermion current and is a scale mass associated with the 

theory (Alexander et al. 2011). In comoving coordinates this equation becomes: 

d\A dadA = j ^-V 2 A- a(t) 2 -^V x A (21) 

v 1 dt 2 dt dt w a(t) M* v 1 

In order to use this in our code, we separated the equations of motion for the gauge field 
and the Friedman equations into a system of first order in time differential equations. 

f - ! < 22 > 

at a 

^ = Ja 3 + V 2 l/a - a 2 -t(V x A) (23) 
dt M* v ' 

d 4 = aH (24) 
dt 

dH 8n . tt9 . . 

IT = T pla ~ H (25) 

The average energy density and pressure are calculated based on equation 15 as 
p = jj J2k=i \-A-k ' Jk\ an d P = —jf J2k=i l^-fc ' <^k\- Here N represents the total number 
of grid points in the computational domain. The scale factor and Hubble parameter 
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therefore depend on the average energy density and not the local field dynamics. Because 
of this, the gravitational back-reaction of the Chern-Simons gauge field on the field 
equations is assumed negligible. 

Once the equations were established we input the initial conditions assuming that 
the current was constant in time (it should be slowly varying) and that the gauge field 
was composed of a random (white noise) spectrum. In order to generate the initial gauge 
field, we used a random number generator to create a random spectrum with amplitude 
up to the calculated maximum amplitude, A , in each direction. The magnitude of the 
gauge field was then held equal to the initial amplitude A . The initial value for the 
other variables is given below. 

\A \ = 10" 5 M P (26) 

J = 10~ 20 Mp(x + y + z) (27) 
f) 

— = 2.18 x l(r 5 Mp (28) 

M* 

a = 1.0 (29) 

H = ^\A -J\ (30) 

— * 

Z Q = H ao x random number(—l, 1) (31) 

The code was then run on the University of Houston Texas Learning & Computation 
Center's Xanadu cluster using a variety of time-steps, grid sizes and resolutions in order 
to obtain consistent results. We ran the code using resolutions ranging from 34 3 to 132 3 
and grid sizes ranging from 10 4 to 10 7 cubed. These grid sizes correspond to the size 
of the horizon throughout the simulation. We adjusted the courant size from a low of 
0.001 to a high of 0.5. These runs utilized anywhere from 16 to 48 processors and the 
runs lasted anywhere from 2 hours to 7 days. We also ran the code using three different 
differencing schemes including 2nd order finite differencing, 4th order finite differencing 
and Fourier spectral differencing. 

Because the initial units were entered as Planck units, we assumed that the physical 
grid (horizon) size corresponded to Planck length and the timing output could be 
interpreted as Planck time. The output was analyzed using several tools including 
ygraph and Pro Fit (www.quansoft.com). We used both Pro Fit and a custom routine, 
which utilized FFTW, to perform a spectral analysis of the gauge field. 



4. Results 



Assuming that the inflationary period begins when the second derivative of the scale 
factor is positive, we find that inflation occurs in our simulation whenever the grid size 
is on the order of 10 6 or larger cubed. In these cases, we find that inflation consistently 
begins around 9.11 x 10~ 38 s. This roughly agrees with the established inflation theory. 
The scale factor continues to increase exponentially until the end of the evolution. 
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Figure 1. The scale factor before and during the inflationary period of the simulation. 



The power spectrum of the gauge field clearly shows the development of an inverse 
energy cascade during the inflationary phase. The initial power spectrum showed a flat 
(random noise) energy distribution with energy variations of up to 2 orders of magnitude. 
As the system evolved, the lowest order modes became more energetic relative to all 
other parts of the spectrum with over 5 orders of magnitude more energy than the highest 
frequency modes. Later in the evolution, the spectrum appears to have flattened out 
as energy returned to the high-frequency modes and variations in the power spectrum 
became smaller. Here high frequency modes correspond to small scale inhomogeneities 
so during inflation, the system appears to become more homogenous in local regions of 
space-time. 

For computational grids with a horizon size of 10 6 planck lengths cubed or larger, we 
obtained results consistent with inflation, regardless of the step-size, differencing method 
or other numerical factors. For smaller grids, however, this exponential expansion does 
not take place. This implies that the inflation process requires wavelengths on the order 
of 10 6 planck lengths. This is most likely a function of our initial conditions for the 
gauge field, current and Chern-Simons interaction terms. 
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Figure 2. The Power Spectral Density of the Gauge Field at different planck times. 
The vertical axis shows the log of the deviation from the mean. The horizontal axis 
shows frequency in planck units. 



5. Discussion and Conclusions 



The simulation presented here clearly supports much of the hypotheses proposed by 
Alexander, Marciano and Spergel: the evolution equations given drive the energy into 
a narrow band of modes and cause the universe to expand exponentially. Though 
the theory does not predict an initial start time for inflation, our value of about 
ti = 10 -37 s is consistent with the time predicted by the most general models of 
inflation (Weinberg 2008). This suggests some universality in the initial conditions 
required for inflation, regardless of the model; in the future, we can study this question 
analytically for this model. However, we do not see inflation end in our simulation, so 
we cannot verify the validity of the mechanism for ending inflation, or predict a stopping 
time numerically. 

The power spectrum in Figure 2 supports the hypothesis that an initially random 
gauge field spectrum becomes concentrated in low frequencies during inflation. This 
analysis allows us to see the progression from an almost completely random spectrum 
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near the beginning to the energy being concentrated in a narrow band of low-frequency 
modes during the inflation event. In addition to the dominance of low-frequency modes, 
the higher modes are suppressed during inflation; the beating pattern probably results 
from imperfections in the random number generator. This reflects the behavior predicted 
in (Alexander et al. 2011). Though we were not able to see the end of inflation in this 
simulation, the latest time power spectrum suggests that inflation is in fact ending. The 
dominance of low-frequency modes has subsided, and we are mostly left with random 
noise behind a flat power spectrum, indicating that the end of inflation is near. 

It is interesting that exponential growth does not occur for horizon sizes that are too 
small. This could be due to the unavailability of very long-wavelength (low frequency) 
modes for the gauge field. In the future, accurate measurements of these frequencies 
could allow us to verify that they are superhorizon modes. However, the most logical 
next step in this research will involve improving the dynamics of the gauge field and 
fermion current to better predict the length and end of inflation. This may involve 
changes to the evolution equations for the gauge field or fermion current during the 
inflation event which may also have an impact on the field dynamics at the sub-horizon 
level. 
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